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' Abstract 



We present results from a simulation of quenched overlap fermions with Liischer-Weisz gauge held 
action on lattices up to 24 3 48 and for pion masses down to ~ 250 MeV. Among the quantities we 
study are the pion, rho and nucleon masses, the light and strange quark masses, and the pion decay 
constant. The renormalization of the scalar and axial vector currents is done nonperturbatively in 
the RI — MOM scheme. The simulations are performed at two different lattice spacings, a ~ 0.1 
fm and ~ 0.15 fm, and on two different physical volumes, to test the scaling properties of our 
action and to study hnite volume effects. We compare our results with the predictions of chiral 
perturbation theory and compute several of its low-energy constants. The pion mass is computed 
in sectors of fixed topology as well. 
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I. INTRODUCTION 



Lattice simulations of QCD at small quark masses require a fermion action with good 
chiral properties. Overlap fermions possess an exact chiral symmetry on the lattice Q], 
and thus are well suited for this task. Furthermore, overlap fermions are automatically 0(a) 
improved if employed properly jjj. 

Previous calculations of hadron observables from quenched overlap fermions have been 
limited to larger quark masses and/or coarser lattices due to the high cost of the simula- 
. To ensure that the correlation functions this involves are not overshadowed 
by the exponential decay of the overlap operator jg|, the lattice spacing a should be small 
enough such that m#a <C 2 for mesons and m#a <C 3 for baryons, where m# is the mass 
of the hadron. In addition, the spatial extent of the lattice L should satisfy L ^> 1/ (2f n ) in 
order to be able to make contact with chiral perturbation theory j^|. 

Over the past years we have done extensive simulations of quenched overlap fermions 
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Furthermore, we have employed overlap fermions to 



probe the topological structure 



13j. In this paper we shall give 



of the QCD vacuum at zero j3| and at finite temperature 
the technical details of our calculations and present results on hadron and quark masses and 
the pseudoscalar decay constant, including nonperturbative renormalization of the scalar, 
pseudoscalar and axial vector currents. The bulk of the simulations are done on the 24 3 48 
lattice at lattice spacing a ~ 0.1 fm. Our results on the spectral properties of the overlap 
operator and nucleon structure functions jn| will be reported elsewhere in detail. 

The paper is organized as follows. In Section II we discuss the action and how it is 
implemented numerically. In Section III we give the parameters of the simulation. In 
Section IV we present our results for the hadron masses and the pseudoscalar decay constant. 
The latter is used to set the scale. We compare our results with the predictions of chiral 
perturbation theory, and attempt to compute some of its low-energy constants. In Section 
V we compute the renormalization constants of the scalar and pseudoscalar density, as well 
as the axial vector current, nonperturbatively, and in Section VI we present our results for 
the light and strange quark masses. Finally, in Section VII we conclude. 
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II. THE ACTION 



The massive overlap operator is defined by 



D=\1-^J D N + m q (1) 
with the Neuberger-Dirac operator Dn given by 

D» = £(l+ , DW{P) ) , D W ( P )=D W -P, (2) 
Dl^Dwip)) a 

where D\y is the massless Wilson-Dirac operator with r = 1, and p G [0,2] is a (nega- 
tive) mass parameter. The operator has n_ + n + exact zero modes, D^ip® = with 
n — 1, • • • , ri_ + ri+ , where n_ (n + ) denotes the number of modes with negative (positive) 
chirality, 75^° = — ip^ (75^° = +V ; n)- The index of D N is thus given by v = n_ — n + . The 
'continuous' modes A.j, D N ipi = A.^, satisfy (^,75^) = and come in complex conjugate 
pairs A-, A*. 

To evaluate it is appropriate to introduce the hermitean Wilson-Dirac operator 
H w (p) = 7 5 Dvf(p) ; such that 

D N = P(l + l5 sgn{H w (p)}) , (3) 
a 

where sgn{if } = H / v ' H 2 . The sign function can be defined by means of the spectral 
decomposition 

sgn{H w (p)} = s g n M XiXl , (4) 

i 

where Xi are the normalized eigenvectors of H\y(p) with eigenvalue /ij. Equation (@J) is, 
however, not suitable for numerical evaluation. We write 

N 

sgn{H w (p)} = Y sgn{/i 4 } Xl xl + P± sgn{H w (p)} , (5) 
i=i 

where 

Pf=l-EXiXl (6) 

i=l 

projects onto the subspace orthogonal to the eigenvectors of the N lowest eigenvalues of 



\H w (p)\, and approximate P^sgn{H w (p)} by a minmax polynomial [JjJ. More precisely, 
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we construct a polynomial P(x), such that 

1 



P(x) 



< e, x E \fii N+1 ,fi 



N+l' FmaxJ ' 



(7) 



where /xjv+i (AWx) is the lowest nonzero (largest) eigenvalue of \P±Hw(p)\- We then have 



N 



sgn{H w (p)} = sgn{^} X t x\ + P?H w {p) P(H 2 w (p)) . 



(8) 



i=i 



The degree of the polynomial will depend on e and on the condition number of H^r(p), 
K = on the subspace {Xi I (1 - P±)Xi = 0}. 

We use the Liischer-Weisz gauge action jl| 



S[U] = 1 



^ - Re Tr (1 - £/ p i aquct te) + Qi - Re Tr (1 - U iec 

parallelogram ) 



^0 

plaquette rectangle 



(9) 



+ c 2 ^ ^ReTr(l-f/ p 

parallelogram 

where t/piaquette is the standard plaquette, {7 re ctangie denotes the closed loop along the links of 
the 1x2 rectangle, and C/paraiieiogram denotes the closed loop along the diagonally opposite 
links of the cubes. The coefficients c\, c 2 are taken from tadpole improved perturbation 
theory 

(10) 



ci (1 + 0.4805a) c 2 0.03325a 
with c + 8ci + 8c 2 = 1, where 



20u 2 



c 



^0 — ( gTr (J/plaquette) 



a 



log(u 



3.06839 



We write 
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co- 



lli) 



(12) 



After having fixed f3, the parameters ci, c 2 are determined. In the classical continuum limit 
u — > 1 the coefficients ci, c 2 assume the tree-level Symanzik values ci = -1/12, c 2 = 0. 
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III. SIMULATION PARAMETERS 



The simulations are done on the following lattices: 



p 


Volume 


r /a 


8.00 


16 3 32 


3.69(4) 


8.45 


16 3 32 


5.29(7) 


8.45 


24 3 48 


5.29(7) 



The scale parameter r /a was taken from [l6j. The couplings have been chosen such that 
the 16 3 32 lattice at (3 = 8.0 and the 24 3 48 lattice at (3 = 8.45 have approximately the same 
physical volume. This allows us to study both scaling violations and finite size effects. 

We have projected out N = 40 lowest lying eigenvectors at /3 — 8.0 and N = 50 (N = 
10) at = 8.45 on the 24 3 48 (16 3 32) lattice. These numbers scale roughly with the 
physical volume of the lattice. The degree of the polynomial P has been adjusted such that 
1/ ^/ H^ v (p) is determined with a relative accuracy of better than 10~ 7 . 

The mass parameter p influences the simulation in two ways. First, it affects the locality 
properties [s| of the Neuberger-Dirac operator. In Fig. ^ we show the effective range of Dn, 




( / max \D N (x, y)\ 



with respect to the Euclidean distance 



\x-y\\=r. 



(14) 



y I u 



\x\ 




(15) 



Asymptotically, F(r) oc exp — fir /a, where p depends (among others) on p. (Numerically, 
p m 2 v, where v refers to the taxi driver distance U) We want p to be as large as possible, 
in particular 2p 3> m#a (3p 3> m#a) for mesons (baryons). Secondly, the condition number 
of PfHyy(p), k = /^max/^TV+i' depends on p as well. In Fig. El we show the p dependence of 
p and k on the 16 3 32 lattice at (3 — 8.45 for N = 10. Test runs show, however, that k does 
not decrease significantly anymore if we increase N further. We have chosen p = 1.4, which 
is a trade-off between a small condition number k and a large value of u. At this value of 
p we find p = 1.11(1), which is consistent with the results obtained in j8j from the Wilson 
gauge action. 
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FIG. 2: Condition number and p, on the 16 3 32 lattice at (3 = 8.45 for p = 1.2, 1.3, 1.4, 1.5 and 
1.6, from left to right. 



6 



The simulations are performed at the following quark masses: 



p 


V 


am q 


8.00 


16 3 32 


0.0168 


0.0280 


0.0420 


0.0560 


0.0840 


0.1400 


0.1960 


8.45 


16 3 32 








0.0280 


0.0560 


0.0980 


0.1400 


8.45 


24 3 48 




0.0112 


0.0196 


0.0280 


0.0560 


0.0980 


0.1400 



(16) 



This covers the range of pseudoscalar masses 250 < mps < 900 MeV as we shall see. The 
lowest quark mass was chosen such that mpsL > 3 (L being the spatial extent of the lattice). 
On all our lattices we have L 3> l/(2f n ). 

0(a) improvement, both for masses and on- and off-shell operator matrix elements, is 
achieved by simply replacing D by 



D 



imp 



am\ ( a 

) D [1 D 

2p) V 2p 



(17) 



in the calculation of the quark propagator. Apart from the multiplicative mass term, this 
amounts to subtracting the contact term from the propagator. In the following we shall 
always use the improved propagator, without mentioning it explicitly. The eigenvalues of 
lie on a circle of radius p/a around (p/a, 0) in the complex plane, while the eigenvalues 

of the improved operator -D™ p = ( 1 — ^~D N ) lie on the imaginary axis. 



2p 



IV. HADRON MASSES AND PSEUDOSCALAR DECAY CONSTANT 



Let us now turn to the calculation of hadron masses and the pseudoscalar decay constant. 
Before we can compare our results with the real world, we have to set the scale. We will use 
the pion decay constant to do so, for reasons which will become clear later. The pion decay 
constant derives from the axial vector current, which has to be renormalized in the process. 



A. Calculational Details 

The coefficients ci, c 2 of the gauge field action are jl(| c\ = —0.169805, c 2 = —0.0163414 
at (3 = 8.0 and c x = -0.154846, c 2 = -0.0134070 at /3 = 8.45. For the gauge field update 
we use a heat bath algorithm, which we repeat 1000 times to generate a new configuration. 
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The inversion of the overlap operator D is done by solving the system of equations 

Ax = y, (18) 

where A = D^D and y is the relevant source vector. We use the conjugate gradient algorithm 
for that. The speed of convergence depends on the condition number of the operator A, 
k(A) = ^max/^minj where z/ max (z/ min ) is the largest (lowest) eigenvalue of A. For reasonable 
values of the quark mass we have k(A) oc 1/m 2 . Thus, the number of iterations, no, needed 
to achieve a certain accuracy will grow like ri£> oc l/m q as the quark mass is decreased. 

The convergence of the algorithm can be accelerated by a preconditioning method. In- 
stead of (|18|) we solve the equivalent system of equations 

ACx = Cy = Ax, (19) 

where C is a nonsingular matrix, which we choose such that k(A) <C k(A). Our choice is 

C=l + £(i-lWj, (2°) 
i=i v * J 

where ivi) are the normalized eigenvectors (eigenvalues) of A. The condition number of the 
operator A is by a factor v n +\lv\ smaller than the condition number of the operator A, and 
the number of iterations in the conjugate gradient algorithm reduces to n D oc 1/ ^Jv n +i + m 'q, 
which depends only weakly on the quark mass m q . We have chosen n = 80, and the inversion 
was stopped when a relative accuracy of 10~ 7 was reached. 

In the calculation of meson and baryon correlation functions we use smeared sources to 
improve the overlap with the ground state, while the sinks are taken to be either smeared 
or local. We use Jacobi smearing for source and sink To set the size of the source, we 
have chosen k s = 0.21 for the smearing hopping parameter and employed N s = 50 smearing 
steps. 

To further improve the signal of the correlation functions, we have deployed low-mode 
1^ | in some cases by breaking the quark propagator into two pieces, 



averaging 



- ^i{x)tp\{y) 



\ im p i m ' ^ ' 



where the sum extends over the eigenmodes of the ne lowest eigenvalues of -D™ p , and the 
remainder. The contribution from the low-lying modes (|2~Tjl is averaged over all positions of 
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the quark sources. As the largest contribution to the correlation functions comes from the 
lower modes, we may expect a significant improvement in the regime of small quark masses. 
We have chosen ri£ = 40, mainly because of memory limitations. 



B. Lattice Results 



The calculations are based on 900 — 1300 gauge field configurations at the lowest four 
quark masses at (3 — 8.0, and on 200 — 300 configurations elsewhere. We consider hadrons 
only with all quarks having degenerate masses. 

Pion Mass 

To compute the pseudoscalar mass, mps, we looked at correlation functions of the pseu- 
doscalar density P = ip^^ip and the time component of the axial vector current A4 = ^7475^ ■ 
In Fig|3]we show the corresponding effective mass for our four lowest quark masses on the 
24 3 48 lattice. Local sinks are found to give slightly smaller error bars than smeared sinks, 
so that we will restrict ourselves to this case. Both correlators give consistent results. We 




t/a 



FIG. 3: The effective pseudoscalar mass from the correlation function of the axial vector current 
A4 on the 24 3 48 lattice at j3 = 8.45, using smeared sources and local sinks. 
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V 


am q 


amps 


amy 




afps 


m PS [MeV] 






0.0168 


0.190(1) 


0.643(5) 


0.793(5) 


0.075(1) 


239(1) 






0.0280 


0.235(1) 


0.64935) 


0.821(4) 


0.076(1) 


295(1) 






0.0420 


0.281(1) 


0.65923) 


0.863(3) 


0.078(1) 


353(1) 


8.00 


16 3 32 


0.0560 


0.321(1) 


0.669(2) 


0.890(3) 


0.080(1) 


403(1) 






0.0840 


0.388(1) 


0.695(3) 


0.952(7) 


0.082(1) 


488(1) 






n i4nn 


0.502(1) 


n 7^1 (o\ 


i r\7d(7\ 

J- .U t t 1 




UOl I ± J 






0.1960 


0.599(1) 


0.815(1) 


1.188(7) 


0.097(1) 


753(1) 






0.0280 


0.212(3) 


0.441(6)* 


0.595(6)* 


0.053(1) 


396(8) 




16 3 32 


0.0560 


0.289(2) 


0.482(4)* 


0.675(4)* 


0.058(1) 


545(4) 


8.45 
















0.0980 


0.384(2) 


0.537(4) 


0.784(7) 


0.064(1) 


727(4) 






n i Ann 


0.467(2) 


0.595(3) 


0.886(6) 


0.070(1) 


883(4) 






0.0112 


0.139(1) 


0.429(6)* 


0.551(12)* 


0.051(1) 


264(4) 






0.0196 


0.177(1) 


0.442(6)* 


0.572(11)* 


0.052(1) 


336(2) 




24 3 48 


0.0280 


0.209(1) 


0.452(3)* 


0.600(10)* 


0.054(1) 


396(2) 


8.45 








0.0560 


0.292(1) 


0.481(3) 


0.674(12) 


0.058(1) 


551(2) 






0.0980 


0.388(1) 


0.538(2) 


0.788(11) 


0.065(1) 


731(2) 






0.1400 


0.412(1) 


0.597(1) 


0.892(11) 


0.071(1) 


887(2) 



TABLE I: Hadron masses and pseudoscalar decay constant. The numbers marked by * are obtained 
with low-mode averaging. To convert mps to physical units, we have used the result in (|28|). The 
error on mps in the last column is purely statistical. 

will use the results from the axial vector current correlator here, because it results in a wider 
plateau as the pseudoscalar correlator, in particular at the larger quark masses. We fit the 
correlator by the function A cosh {mpsit — T/2)), where T is the temporal extent of the 
lattice, over the region of the plateau. The results of the fit are listed in Table 

Rho and Nucleon Mass 

To compute the vector meson mass, my, we explored correlation functions of opera- 
tors Vi = ip^iip and = ^Ci^) {% = 1, 2, 3). We found that the operator Vi, in combination 
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FIG. 4: The effective nucleon mass on the 24 3 48 lattice at (3 = 8.45, using smeared sources and 
local sinks. The horizontal lines indicate the fit interval as well as the value and error of the mass. 
The data points at the lowest three quark masses have been computed with low-mode averaging. 
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FIG. 5: The same as the previous figure, but for the lowest three quark masses without using 
low-mode averaging. 
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with a local sink, gives the best signal. 

For the calculation of the nucleon mass, m^, we used = Sabc^^C^tp ) (where 
C = 7472) as our basic operator, where we have replaced each spinor by xjj — > tp NR = 
(1/2)(1 +74)^ [18]. These so-called nonrelativistic wave functions have a better overlap 
with the ground state than the ordinary, relativistic ones. In Fig. 0] we show the effective 
nucleon mass for all our six quark masses on the 24 3 48 lattice, where for the lowest three 
quark masses we have employed low-mode averaging. We find good to reasonable plateaus 
starting at t/a > 8. In Fig.EJwe show, for comparison, the result obtained without low-mode 
averaging. In this case the situation is less favorable. The nucleon mass is obtained from a 
fit of the data by the correlation function ylexp(— m^t) + B exp(— m^* (T — £)), where m^* 
is the mass of the backward moving baryon, over the region of the plateau. 

The results for the rho and nucleon masses are listed in Table H] Note that amy <C 2/j, 
and arriN <C 3/z, respectively, are satisfied in all cases. In Fig. El we show an APE plot for 
our three lattices. At our smallest quark masses we have mps/my « 0.3. The APE plot 
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1 



m- 
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_i 1 i_ 



_i 1 i_ 



0.2 0.4 0.6 

(m ps /m v ) 2 



0.8 



FIG. 6: APE plot on the 24 3 48 lattice at = 8.45 (O) and on the 16 3 32 lattices at /? = 8.0 (□) 
and (3 = 8.45 (a). 
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shows no scaling violations outside the error bars and no finite size effects. 
Pion decay constant 

The physical pion decay constant is given by 

(0\A 4 \it) = m„U, (22) 

where Ay, is the renormalized axial vector current, A^ = Z^A^. Using the axial Ward 
identity 

d li A„ = 2m q P, (23) 

where P is the local pseudoscalar density, and considering the fact that m g P is a renormal- 
ization group invariant, we obtain 

/. = ^<0|Pk>. (24) 

On the lattice we consider the correlation function 

(P s (t)P s '(0)) = —^(Q\P°\PS)(PS\P s '\0) [exp(-m PS t) + exp(-m PS (T - *))] 

2am PS (25) 

= A ss ' [exp(-m PS t) + exp(-m P5 (T - t))} , 

where the superscripts s, s' distinguish between local (L) and smeared (S) operators. From 
this we obtain 

afps = am \-^rj 7w (26) 

We thus find afps by computing A LS and A ss . In Table H] we give our results. In our 
notation the experimental value of f w is 92.4 MeV. 

Comparing our data on the 16 3 32 and 24 3 48 lattice at (3 — 8.45 in Table H] piece by piece, 
we also find no finite size effects down to the lowest common pseudoscalar mass. 



C. Setting the Scale: Pion Decay Constant 

We will use the pseudoscalar decay constant to set the scale. The reason is that fps is 
an analytic function in mp s for degenerate quark masses j^], in contrast to my and m^, 
which exhibit nonanalytic behavior. We thus expect that f P $ extrapolates smoothly to the 



13 



0.12 



0.1 



0.08 - 



m 

Oh 



05 



0.06 



0.04 



0.02 - 




_i i i_ 



_i i i_ 



0.1 0.2 0.3 

(am ps ) 2 



0.4 



FIG. 7: Chiral extrapolation of the pseudoscalar decay constant on the 24 3 48 lattice at /3 = 8.45 
(O) and on the 16 3 32 lattice at [3 = 8.0 (□). 



chiral limit. In quenched chiral perturbation theory |22j, |23j to NLO we have 1 



f PS = fo[l + al mps 



0(m 



PS) 



(27) 



'2(47r/o) 2 , 

In Fig. [7|we show our data together with a quartic fit in the pseudoscalar mass. The lattice 
spacing is obtained from requiring fpg = fn = 92.4 MeV at the physical pion mass. Using 
the ro/ a values given in ()13|) . we can convert the lattice spacing a into the dimensionful scale 
parameter r . Altogether, we obtain 



(3 


V 


afo 


a t 


a [fm] 


fo [MeV] 


r [fm] 


8.00 


16 3 32 


0.073(1) 


1.5(4) 


0.157(3) 


92(1) 


0.58(2) 


8.45 


24 3 48 


0.049(1) 


1.9(4) 


0.105(2) 


91(2) 


0.56(2) 



(28) 



Note that a\, fo and ro come out independent of the lattice spacing within the error bars, 
which, once more, indicates good scaling properties of our action. The coefficient ot\ turns 



1 Here and in the following we shall adopt the notation a\ = 1287r 2 Lf , L\ being the conventional Gasser- 
Leutwyler coefficients [2^. The superscript q stands for quenched. 
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FIG. 8: Chiral extrapolation of the pseudoscalar mass on the 24 3 48 lattice at = 8.45 (O) and 
on the 16 3 32 lattice at (3 = 8.0 (□). The curves show the fits for a\ = 0. 



out to be in agreement with the phenomenological value of 1.83 (L\ = 0.00145) reported 
in 2^. 



D. Comparison with Chiral Perturbation Theory 



We shall now compare our results for the pseudoscalar, vector meson and nucleon mass 
with the predictions of chiral perturbation theory and attempt to extrapolate the lattice 
numbers to the chiral limit. 

Pion Mass 

We plot the pseudoscalar masses as a function of the quark mass in Fig. |HJ Quenched 



chiral perturbation theory 



m 



PS 



A<1 



22j predicts in the infinite volume 
2 



21 



26] 



S 



-oil 



y (In y + l) + 



(2ag - a\ 



Oil 



(29) 
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with 

y = ^r ( 3 °) 

x 

where A = 2£//q, E being the 'bare quark condensate', and A x denotes the scale at which 
the «j's are being evaluated. The traditional value is A x = 47rfn, which we will also adopt 
here. For the parameter 5 chiral perturbation theory predicts (22, I^J 

x _ $ (31) 



with /1q = m^i + m 2 — 2m 2 K = (870 MeV) 2 . This gives 5 = 0.183. The parameters fo and a\ 
are known from our fit of fps and are given in (|28|). 

A much sought after quantity is the parameter 5. Though unphysical, it would be a great 
success of the calculation, and of quenched chiral perturbation theory as well, if 5 turned 
out to be in agreement with the predicted value. We shall try to determine 5 directly from 
the data. Let us write 



2 2 

m PS „.. m PS 



z=^, w=^ (32) 
A x m q 

and introduce the effective 6 parameter 

! \nz'w-lnzw' 

d cS = 1 + / > I 33 ) 

w — w 

where z, z' and w, u;', respectively, are adjacent data points. It is easy to see that 

lim 5 cS = 5 . (34) 

In Fig. El we show S e f[ as a function of the quark mass. In the case of our high statistics run 
on the 16 3 32 lattice at f3 = 8.0 we are able to extrapolate 5 c q to the chiral limit. We obtain 
5 = 0.18(4), in agreement with the prediction of quenched chiral perturbation theory. On 
the 24 3 48 lattice at (3 = 8.45 our current statistics does not allow such an extrapolation. 
But the data for S e s are not inconsistent with the predicted value of 5. 

The Witten-Veneziano formula relates /ig to the topological susceptibility 

Xt = — > ( 35 ) 
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FIG. 9: The parameter 5 e g on the 16 3 32 lattice at (5 = 8.0 together with a linear fit (top) and on 
the 24 3 48 lattice at (5 = 8.45 (bottom), as a function of the average quark mass m q = (m q + m' q )/2. 
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FIG. 10: The pseudo-scalar mass on the 16 3 32 lattice at (3 = 8.0 for \Q\ = 0, • • • , 9, from left to 
right. The data have been displaced horizontally. The true quark masses are indicated by the 
arrow at the bottom rim of the figure. The curve is from Fig. |SJ 



where Q is the topological charge and V the lattice volume. The result for S is 

1 



Xt , (36) 



which suggests that the pseudoscalar mass depends on the topological charge \Q\. This 
turns out to be indeed the case. In Fig. Elwe show the pseudoscalar mass for various charge 
sectors, where the charge Q is given by the index v of Dm. We observe a strong increase of 
5 with increasing \Q\, and contrary to the findings in [28], we do not expect the effect to go 
away in the limit V — > oo, Xt fixed. It would be interesting to search other quantities for a 
| Q | -dependence as well. 

Let us now turn to the fit of (|29jl to the data. Knowing fo and a|, this leaves us with four 
free parameters. Because our data do not allow an uncorrelated fit of all four parameters, 
we have to make a choice and fix one of them. We consider two cases. In the first case we 
fix a$ at 0, while in the second case we fix 5 at its theoretical value of 0.183. The two fits 
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-0.2(1) 


1.4(4) 


0.8 



(37) 



where we have omitted the heaviest mass point at /3 — 8.00. The numbers shown in italics 
are the numbers that we fixed. It is not expected that A scales. Assuming 5 = 0.183 
and taking /„ from flUJ), we obtain a 3 S = 0.0039(1) at = 8.0 and a 3 S = 0.00138(5) at 
(3 = 8.45, respectively. We shall return to £ and the fit function (J29*|) when we compute the 
renormalized chiral condensate and quark masses. Combining the results on both lattices, 
we obtain ag = 1.5(4) for 5 = 0.183. This is to be compared with 2f]] a s = 0.8(4) in full 
QCD. In Fig. Owe compare the fits with the data. 

Rho Mass 

In Fig. ^2 we plot the vector meson masses as a function of the pseudoscalar mass, where 
we have used the results of (J28|) to convert the lattice numbers to physical values. Quenched 
chiral perturbation theory predicts 



m y = Cq + CX/2 m PS + C\ m 2 p S + C^ /2 m 3 PS 



(38) 



where mps is the lattice pseudoscalar mass as described by (J29|) . The coefficient Cy 2 is 
expected to be negative, so that the chiral limit is approached from below. Our data show 
no indication of a cubic term, and so we shall drop that. A quadratic fit in the pseudoscalar 
mass gives 
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0.48(6) 



(39) 



Our high statistics run at (3 = 8.0 gives indeed a negative value for Cu2, but perhaps of 
lower magnitude than expected [29(, while at f3 = 8.45 our statistics is not high enough to 
make any statement. The fits are shown in Fig. El One might think that at the lighter 
quark masses one is seeing the lowest two-pion state instead of the rho. In Fig. E]we also 
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FIG. 11: Chiral extrapolation of the vector meson mass on the 24 3 48 lattice at j3 = 8.45 (O) and on 
the 16 3 32 lattice at j3 = 8.0 (□), together with the experimental value (*). The solid curves show 
the fits. The dashed curve in the top left corner shows the energy of the state of two pseudoscalar 
mesons. 

show the energy of two pseudoscalar mesons at the lowest nonvanishing lattice momentum 2 , 
\p\ = 2n/(aL), assuming the lattice dispersion relation to hold. We see that the lowest 
two-pion energy lies well above the vector meson mass because of the finite size of our 
lattice. 

Nucleon Mass 

We plot the nucleon masses as a function of the pseudoscalar mass in Fig. [T21 Quenched 
chiral perturbation theory predicts 



m N = C£ + Cf /2 m PS + Cf m' PS + Cfc m% + ■ ■ 



■iN 
'1/2 



N _2 



where 



a 



» 2 = --(3F-D)\6. 



(40) 



(41) 



Note that the pions in the rho are in a relative p wave. 
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FIG. 12: Chiral extrapolation of the nucleoli mass on the 24 3 48 lattice at f3 = 8.45 (O) and on the 
16 3 32 lattice at j3 = 8.0 (□), together with the experimental value (*). The curves show the fits. 

Assuming the tree-level values F = 0.50 and D = 0.76, we expect Cy 2 = —2.58 5. For 
the theoretical value 5 = 0.183 this would give C^ 2 = —0.47. Of course, F and D may 
be different in the quenched theory. In the N c — > oo limit, for example, F/D = 1/3 giving 
C^ 2 = 0. Again, our data show no indication of a cubic term, and we shall drop that here 
as well. A quadratic fit in the pseudoscalar mass gives 
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At (3 = 8.0 we find some evidence for nonanalytic behavior, but with positive coefficient 
Cy 2 . The fits are shown in Fig. 1121 

Both, the vector meson and nucleon masses scale, within the error bars, with the inverse 
lattice spacing set by the pion decay constant f n . 
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V. NONPERTURBATIVE RENORMALIZATION 



We shall now turn to the determination of the renormalization constants Zs, Zp and Za of 
the scalar and pseudoscalar density and the axial vector current, respectively, which we will 
need in order to compute the renormalized quark mass. We shall employ the RI' — MOM 
scheme j^. Our implementation of this method is described in 32^ . 

We consider amputated Green functions, or vertex functions, Tq, with operator insertion 
O = S, P and A± in the Landau gauge. Defining renormalized vertex functions by 



r%(p) = z q (ri- 1 z (jjL) r (p) 



(43) 



where /i is the renormalization scale, we fix the renormalization constants by imposing the 
renormalization condition 



12 



Tr[rg( P )r o ; Born ] 



1 . 



That is, we compute the renormalization constants from 



Z q ^r l Z ^)^T[Y { P )Y ] Bovn ] 



(44) 



(45) 
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FIG. 13: The renormalization constant Za on the 24 3 48 lattice at (3 = 8.45 at the smallest quark 
mass. 
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FIG. 14: Chiral extrapolation of Za on the 24 3 48 lattice at j3 = 8.45, together with its value in 
the chiral limit (#) . 

With T^Born = 1, T^Born = 75 and I^Born = 7475- 

The renormalization constant of the axial vector current can be directly determined from 
the axial Ward identity 



2m q (P(t)P(O)) 
A (6VU(t)P(0)> ■ 

The wave function renormalization constant Z q can thus be obtained from and Za 

Z q (fj) = Z A A A . 



(46) 



(47) 



In Fig. EI] w e plot Z A . We find that the r.h.s. of (|4*U|) is independent of t, except for the 
points close to source and sink, as expected. We extrapolate Z A linearly in am q to the chiral 
limit, as shown in Fig. El The final result is 



p 


V 


Z A 


8.0 


16 3 32 


1.59(1) 


8.45 


16 3 32 


1.42(2) 


8.45 


24 3 48 


1.42(1) 



(48) 
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The corresponding fully tadpole improved (FTI) perturbative numbers [33] are Za = 1.358 
at j3 — 8.0 and Za = 1.303 at (3 — 8.45. They lie 15% — 8% below their nonperturbative 
values. 

Let us now turn to the calculation of Ag(p), A P (p) and A^(p). We denote the expressions 
at finite m q by A(p, m q ). Strictly speaking, Ag(p, m 9 ) and Ap(p, m q ) cannot be extrapolated 
to the chiral limit. Due to the zero modes, both As(p,m q ) an d Ap(p,m q ) diverge oc l/m q . 
This is an artefact of the quenched approximation. On top of that, Ap(p,m q ) receives a 
contribution oc S/(m g p 2 ). This term is due to spontaneous chiral symmetry breaking 
l^ . We thus expect the following dependence on the quark mass: 



32. 



C s (r>) 

A s (p, m q ) = — + C?(p) + Ci(p) am q , (49) 

A P (p, m q ) = + ^ + C 3 P (P) + C[(p) am q , (50) 

(am q ) z am q 

A A (p, m q ) = C£{p) + C£(p) am q , (51) 

neglecting terms of 0(m q ). This behavior is indeed shown by the data. In Fig. ^] we plot 
As(p,m q ), Ap(p,m q ) and Aa(p, m q ) for three different momenta, together with a fit of (|49|) . 
()50j) and (joTjl to the data. We identify A s (p), A P (p) and Aa(p) with cf(p), c^(p) and cf(p), 
respectively, from which we derive 

z M = ¥n z *' ^M = ^4^- (52) 

As(m) A P (/i) 
We expect Zs{^) = Zp(n) due to chiral symmetry. To test this relation, we plot the ratio 
As/Ap in Fig. We find good agreement between Z$ and Zp for all momenta. In the 
following we shall make use of a combined fit of As(p,m q ) and Ap(p,m q ), in which we set 

ci(p) = cap). 



We are finally interested in Zs in the MS scheme at a given scale /i. To convert our 



numbers from the RI' — MOM scheme, which we were working in so far, to the MS scheme, 
we proceed in two steps. In the first step we match to the scale invariant RGI scheme, 



Z^i = A ni'-MOM {fi)Zs{fih (53) 



and in the second step we evolve Zg GI to the targeted scale in the MS scheme, 
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FIG. 15: Chiral extrapolation of As (top), Ap (middle) and (bottom) on the 16 3 32 lattice at 
(3 = 8.45 for some representative momenta p = (711,77,2,713,714) in units of 2tt /aL (711,77,2,713) and 
ir/aL (714). 
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FIG. 16: The ratio As(n)/ Ap(/j,) as a function of fi on the 16 3 32 lattice at j3 = 8.45. 

zr^) = ^^r l z§ Gi . (54) 

The matching coefficients A RI '~ MOM {^) and A MS (fx) are known perturbatively to four 
loops 0] • In Fig. El we show Z§ GI . The result is not quite independent of the scale 
parameter \i as it should, but shows a linear decrease in /x 2 for \i > 2 GeV. We attribute 
this behavior to lattice artefacts of 0(a 2 fi 2 ). Indeed, the slope of Zg GI at our two different 
(3 values scales like a 2 to a good approximation. We thus fit the lattice result by 




ryRGI ,LAT n x r< I \2 

Z s — C + Cx(afjL) 



(55) 



and identify the physical value of Z§ GI with Cq. This finally gives 
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1.18(2) 


8.45 


16 3 32 


1.02(1) 



(56) 



The four-loop value for Af s (/j,) has been given in |36j. At fi = 2 GeV it is Af s (2 GeV) 
0.721(10). The error is a reflection of the error of A MS . The nonperturbative result at 
(3 = 8.45 is in good agreement with Zg TI ' RGI from 
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FIG. 17: The scale invariant renormalization constant Z$ GI on the 16 3 32 lattice at (3 
and (3 = 8.45 (□). 

VI. CHIRAL CONDENSATE AND QUARK MASSES 



8.0 (O) 



Having determined the renormalization constant of the scalar density, we may now com- 
pute the renormalized chiral condensate and light and strange quark masses. 

Let us first consider the chiral condensate (ipip)- Strictly speaking, (ipip) is not defined 
in the quenched theory due to the presence of a logarithmic singularity in the chiral limit. 
Nevertheless, we may identify — (ipip) with £ and assume that £ renormalizes like (the finite 



part of) the scalar density. In the MS scheme at /x = 2 GeV we then have 

(^) MS {2 GeV) = -Zf B {2 GeV) £ . 



(57) 



Taking £ from our second fit in (}37]) . where we have fixed 5 to its theoretical value 0.183, 
this leads to 





V 




2 GeV) 
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16 3 32 


-[324(8) 


MeV] 3 


8.45 


24 3 48 


-[296(11 


)MeV] 3 



(58) 



The lower number at the larger (3 
other quenched lattice calculations 

spectral density Q, [j^ , which we will address in a separate publication 



ue is in reasonable agreement with phenomenology and 
391 1 . A better way to determine E is by means of the 
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Let us now turn to the evaluation of the quark masses. We shall assume (129(1 with a\ = 
as the basic functional form for the relation between the quark masses and the pseudoscalar 
mass: 



m PS = Xm q + Ym q In m q + Zm 
b 



(59) 



For nondegenerate quark masses, m a q and m q , chiral perturbation theory gives the result 



(mfe) 2 = x| q Z q 1 



+ m b n 

+ Z 1 q q 



m a + m" 



m a q In — m In 
— mP 



(60) 



with no new parameter. In fact, (|fi()J) reduces exactly to (j5H|) in the limit m q — > m*. We 
fit ()59j) to our data to determine the coefficients X, Y and Z. The light quark mass, 
rri£ = (m u + m d )/2, is then found from 

— Xrri£ + Fm« In me + , 



(61) 



while we compute the strange quark mass from 



rn 2 K+ + m 2 



A* 11 



Z 



Y 



me + m s 



mi In 7?t.£ — m s In m s 
— m s 



(62) 



The result is 
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(63) 



The renormalized quark masses are given by 



Zm fTT'q j 



(64) 

where Z m = I/Z5. Combining the bare quark masses in (J6l?() with the results for Zs in (f5TI() 



and below, we obtain in the MS scheme at fi = 2 GeV 
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(65) 



These results are in good agreement with other nonperturbative calculations of the quark 



masses in the quenched approximation 



0,0, HQ 
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VII. CONCLUSIONS 



The extrapolation to the chiral limit has been a major challenge in lattice QCD. We have 
shown that with using overlap fermions it is possible to progress to small quark masses. 
Here we have simulated pion masses down to m n « 250 MeV on both of our lattices. We 
have made an attempt to compute the low-energy constants of quenched chiral perturbation 
theory, with some success. Our results turn out to be consistent with the predicted and/or 
phenomenological values. To fully exploit the potential of overlap fermions at small quark 
masses, one will, however, need a statistics of several thousand independent gauge field 
configurations. 

The pion mass was found to depend on the topological charge \Q\ at small quark masses. 
No such behavior was found for the pseudoscalar decay constant, but a similar effect is 
expected to show up in the chiral condensate 

Overlap fermions, in combination with the Liischer-Weisz gauge field action, show good 
scaling properties already at lattice spacing a m 0.15 fm, owing to the fact that they are 
automatically 0(a) improved, on-shell and off-shell. This helps to reduce the large numerical 
overhead in the algorithm. 

The calculations performed in this paper test many of the ingredients needed for a sim- 
ulation of full QCD, and thus provide a lesson for future applications. 
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